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The processes !^e+w^p + e~ and i/e+p ^ n + e provide the dominant mechanisms for heating 
and coohng the material between the protoneutron star and the staUed shock in a core-coUapse 
supernova. Observations suggest that some neutron stars are born with magnetic fields of at least 
~ 10^^ G while theoretical considerations give an upper limit of ~ 10^* G for the protoneutron star 
magnetic fields. We calculate the rates for the above neutrino processes in strong magnetic fields of 
~ 10^^ G. We find that the main effect of such magnetic fields is to change the equations of state 
through the phase space of e~ and e"*", which differs from the classical case due to quantization of 
the motion of e~ and perpendicular to the magnetic field. As a result, the cooling rate can be 
greatly reduced by magnetic fields of ~ 10^^ G for typical conditions below the stalled shock and 
a nonuniform protoneutron star magnetic field (e.g., a dipole field) can introduce a large angular 
dependence of the cooling rate. In addition, strong magnetic fields always lead to an angle-dependent 
heating rate by polarizing the spin of n and p. The implications of our results for the neutrino-driven 
supernova mechanism are discussed. 

PACS numbers: 25.30.Pt, 26.50.+X, 97.60.Bw 



In this paper we study neutrino processes in strong magnetic fields and their implications for supernova dynamics. 
Although the detailed mechanism by which massive stars produce supernova explosions is still elusive (see Ref. UJ for 
a review), intense research in the past few decades has led to the following prevalent paradigm. At the exhaustion of 
nuclear fuels, the Fe core of a massive star collapses. When nuclear density is reached, the inner core bounces and a 
shock is launched. As the shock propagates outward, it loses energy by dissociating the Fe nuclei falling through it. 
Eventually, the shock is stalled before exiting the outer core. Meanwhile, the inner core is settling into a protoneutron 
star by emitting v^, Pe, t^p,, i^p,, i^t, and Pr- These neutrinos can exchange energy with the material below the stalled 
shock. The dominant energy-exchange processes are 



The forward processes in Eqs. and ^ heat the material through absorption of Ve and De while the reverse processes 
cool the material through capture of e~ and e"*". The competition between heating and cooling of the material by 
these processes is expected to result in net energy gain for the shock, which then propagates outward again to make 
a supernova explosion. This is the neutrino-driven supernova mechanism 

Unfortunately, the current consensus is that the neutrino-driven supernova mechanism does not work in spherically 
symmetric models One group has shown that this mechanism works in three-dimensional models where spherical 

symmetry is broken by convection (Ej. However, this success has not been confirmed by other groups yet. On the 
other hand, magnetic fields may be generated during the formation of protoneutron stars and in turn affect supernova 
dynamics. Magnetic fields of ^ 10^^ G are commonly inferred for pulsars. Observations also suggest that a number of 
neutron stars, the so-called magnetars, have magnetic fields of ~ 10^^ G (see e.g., Refs. [^J3,0). A theoretical upper 
limit of ~ 10^* G may be estimated for the magnetic fields of new-born neutron stars While strong magnetic 
fields may induce supernova explosions directly through dynamic effects such as jet production [l^, and therefore, 
render the neutrino-driven mechanism irrelevant, the details of this magnetohydrodynamic mechanism have not been 
worked out or understood yet. In this paper we address the effects of strong magnetic fields on supernova dynamics 
still in the context of the neutrino-driven mechanism. In particular, we focus on how such fields affect the microscopic 
processes of heating and cooling the material below the stalled shock. 
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I. INTRODUCTION 
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FIG. 1: A sketch of the region of interest. The neutrinosphere of radius Rv = 50 km effectively defines the surface of the 
protoneutron star. The stalled shock is at an average radius Rs = 200 km. At the gain radius Rg, the rate for heating by 
absorption of Ve and Ve equals that for cooling by capture of e~ and e"*". Heating dominates cooling above Rg. 



We describe the neutrino-driven supernova mechanism without magnetic fields in some detail in Sec. The rates 
of the processes in Eqs. ^ and Q in strong magnetic fields are calculated in Sec. lIIII The implications of these rates 
for supernova dynamics are discussed in Sec. II VI and conclusions given in Sec. IVl 



II. THE NEUTRINO-DRIVEN SUPERNOVA MECHANISM 



In this section we give a more quantitative description of the neutrino-driven supernova mechanism in the absence 
of magnetic fields. As this mechanism has not been fully established yet, we will use parameters typical of current 
models to illustrate the essence of these models rather than focus on the numerical details of a specific model. We are 
interested in times of < 0.5 s after the core collapse, which correspond to the critical period for the neutrino-driven 
supernova mechanism. The region of interest is above the protoneutron star but below the stalled shock as illustrated 
in Fig. n We consider that neutrinos are emitted from a neutrinosphere of radius = 50 km that effectively 
defines the surface of the protoneutron star (see e.g., Fig. 11.1 in Ref. ^3). The stalled shock is taken to be at 
an average radius Rs = 200 km. The material in the region of interest has typical entropies of 10 (in units of 
Boltzmann constant per nucleon) and typical temperatures of several MeV. For these conditions, the material can be 
characterized as a gas of n, p, e^, and 7 (photons). The predominant cooling processes are the reverse reactions 
in Eqs. ^ and |(2J) while the predominant heating processes are the corresponding forward reactions [note that similar 
charged-current processes involving f^t(r) and ^^^(t) are energetically forbidden for the neutrino energies and material 
conditions available in supernovae]. Cooling dominates heating near the neutrinosphere. However, the cooling rate 
decreases much more steeply with increasing radius than the heating rate. These two rates becomes equal at the 
gain radius above which heating dominates. Thus, the heating and cooling rates are crucial to a quantitative 
discussion of the neutrino-driven supernova mechanism. These rates are calculated below. 

We start with a description of neutrino emission by the protoneutron star. The luminosity L^^ is taken to be the 
same as the De luminosity ip^. We assume L^^ — Li^^ = 4 x 10^^ erg s~^ during the epoch relevant for shock revival 
by neutrino heating. As Ve and De are roughly in thermal equilibrium with the matter at the neutrinosphere, their 
luminosities approximately correspond to those of black-body radiation for a Fermi-Dirac neutrino energy distribution 
with zero chemical potential: 



L..^L,^-—[T{R,)]^Rl^ixlO' 

In the above equation, T{Ri,) is the temperature at the neutrinosphere. Throughout this paper, we adopt units where 
the Planck constant h, the speed of light c, and the Boltzmann constant k are set to unity. 

Due to the difference in the interaction of i^e and v,, with the protoneutron star matter, their emission is more 
complicated than implied by the crude estimate in Eq. As there are fewer protons to absorb De than neutrons 
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FIG. 2: A sketch of the geometry for calculating the differential neutrino number density d^riiy/dEi^dQi, at radius r. For a specific 
radial direction f , the differential solid angle dS7^ is defined by the polar angle B^, between f and the direction of the neutrino 
momentum pv (the corresponding azimuthal angle 0^ is not shown) . Only neutrinos emitted with \/l ~ {Rv/r)'^ < cos 9v < 1 
can contribute to d?n^ /AE^AQ,^. 



to absorb v,, , Uf. decouple from the protoneutron star matter at higher temperature and density than u,, . This results 
in a higher average energy for (see below). However, due to the steep temperature and density gradients near the 
protoneutron star surface, the radii for and decoupling are essentially the same. Detailed neutrino transport 
calculations show that Ve and Ve can be considered as having the same luminosity and neutrinosphere but significantly 
different average energies (see e.g., Ref. ^.121). The normalized Vf. and energy distributions can be described by 
functions of the form 

1 

^"^^"^ ^ F2(77.)r3exp[(i?,/T,)-ryJ + l' ^"^^ 

where the subscript v refers to z/e or i/g, Ei, is the neutrino energy, and t;,^ are two positive parameters, and F2{rj^) 
is a specific case of the general Fermi integral F„ {rf) defined as 



FJn)= / ^ dx. (5) 

"^'^ io exp(a;-r;) + l ^ ' 

The parameter is related to the average neutrino energy {E^) as 

. fMr. (0, 

We take -q^^ = 77^^ = 3, {E^^J) = 11 MeV, and (SpJ = 16 MeV. For these parameters, T^^ = 2.75 MeV and Tp^ = 4 
MeV, which are close to T{R^) estimated from Eq. JSJ. 

To calculate the rate of heating by neutrino absorption processes, we need the differential neutrino number density 
per unit energy interval and per unit solid angle d^n^/dE^dfli, at radius r > R^, as measured from the center of the 
protoneutron star. For a specific radial direction f , the differential solid angle dVl^, is defined by the polar angle 0^, 
between f and the direction of the neutrino momentum pj, (see Fig.|21l and by the corresponding azimuthal angle (j)v 
We assume that only neutrinos emitted with -^/l — [R^/rY < cosOy < 1 can contribute to d^n^/dE^dQ^ (see Fig.|21l. 
Thus, 



^ / L,UiE,)/{47:^Rl{E,)), for - [Rjrf < cosO, < 1, 



dE^di^i, \ 0, otherwise. 
The corresponding neutrino occupation number f^{E^) is 



Note that 
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FIG. 3: The heating rate per nucleon (solid curve) and the cooling rate per nucleon qi [dashed {S — 10) and dot-dashed 
(5 = 20) curves] as functions of radius r. The gain radius is at Rg — 137 and 68 km for entropies of S = 10 and 20, respectively. 



where d'^nj^/dE^dfli, is evaluated at r = R^. 

To zeroth order in mjj^ where m^y is the nucleon mass, the cross sections for the forward processes in Eqs. ^ and 
@ are 



— (/ +3g )PeEe, 



(10) 



where the subscript i/N refers to i^e absorption on n or i/^ absorption onp, Gp = (292.8 GeV) ^ is the Fermi constant, 
9c is the Cabbibo angle {cos^ 9c = 0.95), f — 1 and g = 1.26 are the vector and axial-vector coupling coefhcients, 
respectivley, of the weak interaction, and Pe and = \/pI + ml are the momentum and energy, respectively, of the 
electron or positron in the final state (we is the electron mass) . To the same order in m^^ , conservation of momentum 
and energy gives 



E^±A, 



(11) 



where A = 1.293 MeV is the neutron- proton mass difference, the plus sign is for absorption of Vf,, and the minus sign 
is for absorption of C'e. Note that there is a threshold energy of A -f- rrie for i^g absorption on p. 
At radius r > R^, the heating rate per nucleon cf^^ due to absorption of v^, and v,, is 

d^n,. 



.(0) 

1h 



oo 

(0) 



dE,, 



■dE. 



lA+m 

(3.83 X lO^Yn + 5.32 x 10%) 



dEi^^drii.^ 

1 - Vi - {R./ry 



MeV s ^ nucleon 



(12) 



where the numerical coefficients correspond to the parameters Lj,_,, Lp^, T^^, Tp^, rj^,^, and rjp^ adopted above, and Yn 
and Yp are the neutron and proton number fractions, respectively. As the material is neutral, Yn ^ l — Y^ and Yp = Yg, 
where Yg is the net electron number per nucleon. The radial dependence in Eq. (|12|l comes from the integration over 

dfip and accounts for the geometric dilution of the neutrino number density. Note that for r ^ R^, g^*^^ cx r~^. 

Taking Ye — 0.5, we calculate q^^'' from Eq. (|12|) and show the result as the solid curve in Fig.|3| 

Next we calculate the cooling rate due to the reverse processes in Eqs. ^ and It is convenient to define a 
volume reaction rate, which gives the rate of e.g., e"*" capture per neutron when multiplied by the e"*" number density 
7ie+ . In the absence of magnetic field, the volume reaction rate V^^)j is simply va^^^ , where v is the relative velocity 
between the reacting particles and af^ is the cross section for e+ capture on n or e~ capture on p. To zeroth order 
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^, the volume reaction rates for these processes are 
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In Eq. lO, 

E,=Ee± A, (14) 

where the plus sign is for e"*" capture on n and the minus sign is for e~ capture on p. The coohng rate per nucleon 
qi^^ is then 



7r2 7o exp[(£;e+/T) + Tye] + 1 TT^ J^A^-ml exp[(£'e-/r) 77e] + 1 



where T is the temperature, yye is the electron degeneracy parameter, and — is the threshold momentum for 

e~ capture on p. 

To evaluate qi^^ we need T and rj^. We take 

, , R,, ( 50 km\ 
T(r) = T{R^)— = 4 MeV (16) 



r 



for Ri, < r < Rs (see e.g., Ref. pj)- We assume that the material in the region of interest can be characterized by a 
typical electron fraction Yg and a typical entropy per nucleon S. We then obtain rj^ together with the matter density 
p from the equations of state: 

= n^- - 71e+, (17) 

mN 

S — Sn + 'S'7 + Sg- + Sg+, (18) 

where Sn, Sj, S^- , and 3^+ are the contributions to S from nucleons, photons, electrons, and positrons, respectively. 
The expressions for n^-, ng+, Sn, S^, S^-, and 3^+ are given in Appendix IXI We note that for extremely relativistic 
e~ and e+, 

ne--n.. = + (19) 



While we always use the expressions in Appendix^to calculate the results presented in this paper, we find that Eqs. 
((T^ and are excellent approximations to the corresponding general expressions (without magnetic fields) for the 
conditions in the region of interest. 

Numerical models 3, .4] show that and S tend to rise sharply over a short distance above the neutrinosphere and 
then stay approximately constant. As cooling always dominates near the neutrinosphere, the gain radius lies in the 
region where Yg and S can be taken as constant. Thus, we assume no radial dependence for Yg and S in determining 
the gain radius. Taking Y^ = 0.5 and = 10 and 20, we calculate p and ije as functions of r from Eq. (|16|l and the 
equations of state [see Eqs. (|17|) and (|18|l and Appendix IX| . The results are shown as the solid {S = 10) and dashed 

{S = 20) curves in Fig.0^ for p and Fig. 0)3 for r]e. We also calculate the corresponding gc"'' from Eq. H15|l and show 
the results as the dashed {S = 10) and dot-dashed {S — 20) curves in Fig. |21 along with the result for q^^^ (solid 
curve). It can be seen that the gain radius is at Rg = 137 and 68 km for 5 = 10 and 20, respectively. As the shock is 
at radius Rs ~ 200 km, there is a large region for net heating below the shock in both cases. Note that the location 
of the gain radius is rather sensitive to 5. A gain radius below the shock {Rg < Rg) exists only for 5 > 10. We have 
also done calculations for different values of Yg and found that the effects of strong magnetic fields to be discussed 
are qualitatively the same for Y^ > 0.3. For clarity of presentation, we focus on the results for Y^ = 0.5. 

In the above calculation of q^'^'' and qi^\ we have ignored Pauli blocking of the final states by the e^, e+, z^e, and z/g 

(see e.g., Ref. 14]) in the region above the neutrinosphere. In the case of ql^\ the e~ and e+ produced by v,, and z/g 
absorption have typical energies of ^ 20 MeV that are much higher than the average energies of and e+ in the gas. 
So Pauli blocking is unimportant for calculating qi^\ especially in the region near and above the gain radius. In the 
case of qc^\ the neutrino occupation number fu{E^) can be calculated from Eq. (jHJ. For the adopted parameters, we 
find fuAEuJ < 0.37 and fuA^uJ < 0.083 for ^1 - {R^/ry < cos 6^ < 1 and UA^uJ = fuA^uJ = otherwise. As 
the range of cos^^ for finite fi,{Ei,) diminishes with increasing r, Pauli blocking in this case is also insignificant in the 
region near and above the gain radius. Thus, ignoring Pauli blocking in the calculation of q^i^^ and qi^"^ has little effect 
on our discussion of the gain radius above. For the same reasons, we will ignore Pauli blocking in the calculation of 
the heating and cooling rates in strong magnetic fields. This approximation will not affect the comparison of the gain 
radii for the cases without and with strong magnetic fields. 
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FIG. 4: The matter density p (a) and the electron degeneracy parameter r]e (b) as functions of radius r. The solid {S — 10) 
and dashed (5 — 20) curves are for the case of no magnetic field. The dot-dashed {S = 10) and dot-dot-dashed (5 = 20) curves 
are for the case of a uniform mae netic field of B = lO^'' G. The magnetic field greatly reduces rje for an entropy of 5 = 10. For 
S — 20, rje is very small even for the case of no magnetic field. 



III. HEATING AND COOLING RATES IN MAGNETIC FIELDS 



In this section we calculate the rates of heating and cooUng by the processes in Eqs. and (0) in strong magnetic 
fields. We consider a uniform magnetic field of constant strength B in the z-direction. Observations indicate that 
neutron stars may have magnetic fields up to ~ 10^^ G long after their birth in supernovae. This suggests that 
magnetic fields of at least ~ 10^^ G can be generated during the formation of some neutron stars. An upper limit of 
^ 10^* G for protoneutron star magnetic fields can be estimated by equating the magnetic energy to the gravitational 
binding energy of a neutron star In this paper we consider protoneutron star magnetic fields of i? ~ 10"'^^ G. 

An obvious effect of the magnetic field is polarization of the spin of a nonrelativistic nucleon due to the interaction 
Hamiltonian 

-fi B = -fia^B. (21) 

In Eq. (|21ll . fi = /xcr is the nucleon magnetic moment, where /x = 2.79^n for p and /i ~ — 1.91/iAr for n with 
fJ-N ^ e/{'2mp) being the nuclear magneton, and cr refers to the Pauli spin matrices. For a nucleon gas of temperature 
T, the net polarization x is 



^ ^ eM^^B/T) ~ exp(-Mi?/r) ^^2) 



cxp{fiB/T) + exp{-fiB/T) 



For \fJ.B/T\ < 1, Eq. ^ reduces to 



3^15 X 10- ,23, 

In addition, the motion of a proton in the xy-plane perpendicular to the magnetic field is quantized into Landau 
levels (see e.g. Ref. J^) with energies 

i?(n.,M = ^+(". + 0^,n. = 0,1,2,..., (24) 

where kp^ is the proton momentum in the z-direction. As eB/nip = 63(i?/10^^ G) keV, a proton in a gas of 
temperature T > 1 MeV can occupy Landau levels with 3> 1 for B 10^^ G. From the correspondence principle, 
the proton motion in this case can be considered as classical. Thus, we only need to take into account polarization of 
the spin by the magnetic field for both p and n. 

For the conditions of interest here, and e+ are relativistic. Their Landau levels |l6| have energies 

Ee{n, pe.) = VpL + 2neB, n = 0, 1, 2, . . . , (25) 
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where Pez is the momentum of e~ or e+ in the z-direction. The result in Eq. H25(l takes spin into account. Note that 
the Landau levels of e~ and e"*" have degeneracy (7„ = 1 (corresponding to a single spin state) for the n = state but 
Qn = 2 (corresponding to two spin states) for all n > states. For a given E,,^ the maximum value rimax of n is 



2eB 



int 



^.45 X 10 



-3 



Ef. 



MeV" 



10^^ G 
B 



(26) 



int 



where [ ]int denotes the integer part of the argument. Thus, and e"*" with Ee < 10 MeV can only occupy the 
Landau level with n = in magnetic fields of B ^ 10^^ G. The phase space of e~ and e+ in this case is dramatically 
affected by magnetic fields. In general, the integration over the phase space is changed from the classical case to the 
case of Landau levels according to 



(27) 



where 51 is the solid angle in the classical momentum space and Pez is restricted to positive values. 



A. Heating Rate in Magnetic Fields 



We first calculate the heating rate due to the forward processes in Eqs. ^ and Q in magnetic fields. As we only 
consider magnetic fields of i? 10^^ G, the energy scale ^/eB = 7.69(i3/10^^ G^^^ MeV is much lower than the mass 
of the W (80 GeV) or Z boson (91 GeV). Thus, the weak interaction is unaffected by such magnetic fields. However, 
such strong magnetic fields can change the cross sections for the forward processes in Eqs. Q and ^ by polarizing 
the spin of n and p in the initial state and by changing the phase space of the e~ and in the final state. The new 
cross sections to zeroth order in mjj^ are derived in Appendix IbI using the Landau wavefunctions of e~ and e+. The 
results are 



if±g)g 
^XTTTl-^cose^ 



.(2) 



r f2 



P + 



cos 9^ + 2x 



if T 9)9 



where we have factored out two energy-dependent terms 
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(28) 

(29) 
(30) 



In Eq. H28|l . 0^ is the angle between the neutrino momentum and the magnetic field, the upper sign is for absoprtion 
on n, and the lower sign is for absoprtion on p. As in the case of no magnetic field, Ee = E^, ± A. 

The angular dependence in Eq. (|28|l is due to parity violation of the weak interaction. The dominant angular 
dependence is associated with polarization of the spin of the n or p in the initial state. The single spin state 
corresponding to the n = Landau level of the or in the final state introduces additional (typically small) 
angular dependence. For comparison, if there is no magnetic field but the spin of the n or p in the initial state is 
polarized, the cross sections for the forward processes in Eqs. and Q are 



Ovn{B = 0) = cr. 



(0) 
vN 



1 



2x^(^cose. 
/ + ^9 



(31) 



In Eq. (|^ . a^^lj = {Gj^ cos^ 9c /t^){P + Sg'^)peEe are the appropriate cross sections for S = and x = as given 
in Eq. (|10|1 . the upper sign is for absoprtion on n, and the lower sign is for absoprtion on p. Note that the 
angular dependence for absorption on p is much weaker than that for v^. absorption on n due to the close numerical 
values of / and g. For illustration, we take x = —0.05 {fj, = -l.QlfiN for n), B = 10^^ G, and calculate a^^n{B) and 
cTu^niB = 0) as functions of E„^. The results are shown as the dotted [(Jv^n{B)] and dashed [tT^^„(i? = 0)] curves in 
Figs. |SJi-c for cos8jy_j — —1, 0, and 1, respectively. For x = 0.05 (// = 2.79/UAr for p) and B = 10^^ G, the angular 
dependence of aj^^p{B) and ay^p{B — 0) is insignificant. The results for cos6p^ = are shown as the dotted [crp^p(i?)] 
and dashed [(Jv^piB — 0)] curves in Fig.|3li. 
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FIG. 5: The cross sections (7^^„(-B) in a magnetic field of B = 10^*^ G (dotted curve) and a^^n{B = 0) (dashed curve) as 
functions of Ve energy E^^ for different angles Qv^ between the directions of the magnetic field and the momentum: (a) 
cosB;^^ = —1, (b) cos 0j/^ = 0, and (c) cos0^^ = 1. A neutron polarization of x = —0.05 is assumed for all cases. The 
cross section (Jv^n{B) smoothed with a Gaussian window function exp [-{E- E^^f/{Q.^ MeVf] (solid curve) oscillates rather 
symmetrically around Gi,^n{B = 0) (dashed curve) at Ei,^ > 20 MeV. (d) The cross sections Gv^p{B) in a magnetic field of 
B = 10'"^ G (dotted curve) and oa^pi^B — 0) (dashed curve) as functions of Ve energy Ej^^ for cos = 0. The dependence 
on cosGi/e is very weak for the assumed proton polarization of x ~ 0.05. The cross section Uv^p{B) smoothed with the same 
window function (solid curve) exhibits the same behavior as in the case of (Ty^„(i3). 



The cross sections (T„^n{B) and ajj^p{B) shown as the dotted curves in Figs. |5ji-d have spikes superposed on a 
smooth general trend. The varying heights of these spikes are artifacts of the plotting tool: all the spikes should have 
been infinitely high as they correspond to "resonances" at — \/ m\ + 2neB, for which a new Landau level opens up. 
These formal infinities disappear when nucleon motion is taken into account [l7| . In practice, these formal infinities 
are effectively smoothed out when integrated over the neutrino energy spectra. To see the behavior of a„N{B) as 
functions of more clearly, we smooth (j^m{B) with a Gaussian window function exp[— (£' — i?,y)^/(0.5 MeV)^]. The 
results are shown as the solid curves in Figs.|3^-d. It can be seen that the smoothed (Ji,^n{B) and ap^p{B) oscillate 
rather symmetrically around the corresponding results for B ^ a.t E^, > 20 MeV. This is because for these neutrino 
energies, the e~ or e+ in the final state can occupy Landau levels with n up to Umax ^ 3 for B — 10^^ G [see Eq. (|26(l ]. 

For rimax > 1, 

^^E^W^=r^^/ ' o ^ <^{neB) ^ E^^W^e, (32) 

,7^0 V - m-e - 2neB Jo yj Ej - - 2neB 

and in the dominant term of a^j\j{B) approaches cr^f^. 

To calculate the heating rate per nucleon qh{B) in magnetic fields, we replace (t^']| with a^NiB) in Eq. Hl'2(l . As 
noted above, the resonances in a^NlB) are smoothed out when integrated over the neutrino energy spectra. This 
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FIG. 6: The ratios {EeO^g'^) / {E^afl) (a) and {E^a^g'^) / {Eea''^}^) (b) as functions of B. The sohd curve is for absorption on 
n and the dashed curve for Ve absorption on p. 



integration results in 



which can be compared with 



E,al^^U{E,)dE,. 



(33) 
(34) 

(35) 



The ratios {E^crg '') / (Eecr^^^) are shown as functions of B for i/g absorption on n (sohd curve) and Ve absorption on p 



(dashed curve) in Fig.|H|i. The corresponding results for {Ee<J g ) / {Eecr^/j) are shown in Fig.ElD. It can be seen that 

for B < 10^^ G, (Eeag^) / (E^al^^) stays close to unity. This is because the dominant contributions to the relevant 
integrals come from with iJ^^ ~ 20 MeV or with E^^ ~ 25 MeV and these neutrino energies correspond to 

n-max ^ 1- The ratio {Egu'^g') / {Eeof^j) is negligible for B < 10^^ G. In general, the overall term involving cr^'' in 

<^un{B) is much smaller than the one involving cr^' [see Eq. (I28|l ]. 

It can be seen from Figs. Et- and 6b that substantial changes to the magnitude of the heating rate only occur for 
B ~ 10^^ G. However, a qualitatively new and quantitatively significant effect already occurs for B ~ 10^^ G. Due to 
polarization of the spin of n or p in the initial state of the heating reactions, the heating rate at position r depends on 
the angle 9 between r and the magnetic field (in the z-direction) . This angular dependence enters through integration 
over the neutrino solid angle (see Fig. [Tj) 



cos 6^ dfijy = ■n[Rylrf' cos B. 



(36) 



B. Cooling Rate in Magnetic Fields 

Next we calculate the cooling rate due to the reverse processes in Eqs. and Q in magnetic fields. The differential 
volume reaction rates dreAr(-B)/dcos for these processes are derived in Appendix IbI The results are 



drejv(B) _ r 

dcos 9i/ 



(0) 



(/±.g)ff 
P + 3g2 



cos 0^ 



P + 3.9- 



■ cos8y + 2x 



(/ T g)g 
P + 3.g2 



(37) 
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FIG. 7: A sketch of the geometry for integration over the neutrino sohd angle to obtain the heating rate per nucleon qh{B) 
in magnetic field. The magnetic field is uniform and in the ^-direction. The cross sections used to calculate qh{B) depend 
on the angle between the directions of the magnetic field (z) and the neutrino momentum (pi/). The integration over the 
neutrino solid angle can be performed by expressing cosOi^ = z • p^^ in terms of 6, 6^, and (j>^ (the azimuthal neutrino angle (j>i' 
corresponding to 6^ is not shown). 



if the e or e+ in the initial state is in the n = Landau level, and 



(0) 



dT,NiB) ^ r 
dcos 9i/ 2 



eN 



(38) 



if the e or in the initial state is in the n > Landau level. In Eqs. ^ and T^^ = [Gj. cos^ ec/{27T)]{p + 
3g^)E^ with i?^ = E'e ± A is the volume reaction rate in the absence of magnetic field as given in Eq. (|13|l . the upper 
sign is for e"*" capture on n, and the lower sign is for e~ capture on p. The dependence of dreAr(i?)/dcos 8,y on the 
direction of the neutrino emitted in the final state again manifests parity violation of the weak interaction. As we 
are not interested in the neutrinos emitted by the cooling processes, we integrate over cos Qi, to obtain the volume 
reaction rate 

TMB) = I + '^(^ ^ '^'/^f' + " = (39) 

[ Tew, n>0. 

Clearly, the volume reaction rates of the cooling processes are not much affected by the magnetic field for |x| <C 1. 

However, magnetic fields also affect the cooling rate through the equations of state for e~ and e"*" (see Appendix 
For a given set of Y^, S, and T, the density p and the electron degeneracy parameter rje in the presence of magnetic 
field differ from those in the case of no magnetic field. Taking — 0.5 and T{r) in Eq. H16|l . we show p{r) for 
B = 10^^ G as the dot-dashed {S = 10) and dot-dot-dashed {S = 20) curves along with the corresponding results for 
B = [solid {S = 10) and dashed {S = 20) curves] in Fig. The comparison for r]e{r) is shown in Fig.^js. It can 
be seen that for the same Ye and T, magnetic fields of i? ~ 10^^ G change (mostly decrease) p slightly but decrease 
r]e greatly for 5* = 10. The same magnetic fields significantly increase p but decrease r]e for = 20. Note that rje is 
already small for B = and S" = 20. 

The cooling rate per nucleon in magnetic fields is 



+Yp7r^y^gn / 77^ J7ff r— j-dpe-z- (40) 
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FIG. 8: Contours of constant qc{B)/qc in increment of 0.1 as functions of S and T for B = 10^® G (a) and as functions of 
B and T for 5* = 10 (b). Magnetic fields of B ~ 10^^ G greatly reduce the cooling rate per nucleon for entropies of 5* ~ 10, 
especially at temperatures of T < 2 MeV. 



In Eq. H40() . the threshold momentum Pe-z,n for e capture on p corresponds to Pg-^„ = — ml — 2neB or 0, 

whichever is larger. Taking Ye = 0.5, we plot contours of constant qc{B)/q"c'^ as functions of S and T for B = 10^® 
G in Fig. and as functions of B and T for S* = 10 in Fig. ^p. It can be seen that magnetic fields of i? ^ 10^^ G 
can decrease the cooling rate significantly for S ^ 10. This is mostly due to the reduction of rje through the effects of 
magnetic fields on the equations of state for and e+ (see Fig. 



IV. IMPLICATIONS FOR SUPERNOVA DYNAMICS 



Now we consider the effects of strong magnetic fields on supernova dynamics using the heating and cooling rates 
discussed in Sec. IIIII As mentioned earlier, observations suggest that magnetic fields of at least ~ 10^^ G can 
be generated during the formation of protoneutron stars. However, little is known about the actual strength and 
topology of protoneutron star magnetic fields. To illustrate the potential effects of such fields on supernova dynamics, 
we consider simple cases of uniform and dipole fields of ^ 10^^ G. 



A. Uniform Field 



We first consider the case of a uniform magnetic field in the z-direction. The heating and cooling rates in such a 
field have been discussed in Sec. IIIII For illustration, we take B = 10^^ G. The corresponding heating rate is 



qh[B ^ 10^6 G) = {3.77 X 10^r„[l - 1.46 x lO^^Xn + (0- 



1.31 X 10"^)$(r) cos 61] 



-5.29 X 10%[1 + 7.45 x lO^^^p - (0.114xp + 7.69 x 10"3)$(r) cos 61]} 

1- Vi-(^.A)' 



MeV s ^ nucleon , 



where Xn and Xp are the net polarization of n and p, respectively, as given in Eq. H22|l . and 



$(r) 



1 - v/1 - {Rjrf 



(41) 



(42) 



Note that <I>(r) increases from 1/2 to 1 as r increases from to r ^ R^. Note also that Xn and Xp are functions 

of r through their dependence on T{r) [see Eqs. (O and For T ~ 2 MeV, Xn 0.03 and Xp 0.04. The 

magnitudes of Xn and Xp increase for lower T. Thus, in the region of T < 2 MeV [r > 100 km, see Eq. (|16ll ]. the 
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FIG. 9: The angle- averaged heating rate per nucleon (solid curve) and the cooling rate per nucleon [dashed {S = 10) and 
dot-dashed (S = 20) curves] in a uniform magnetic field of B = 10^® G as functions of radius r. Compared with Fig.|21 the gain 
radius decreases significantly from 137 km for _B = to 99 km for B = 10^^ G in the case of S = 10 but essentially remains at 
68 km in the case of 5 = 20. 

heating rate in Eq. (|41() varies by at least several percent over — 1 < cos 9 < 1 for a given r. This variation can induce 
or amplify anisotropy in the bulk motion of the material below the stalled shock, eventually producing an asymmetric 
explosion. The protoneutron star would then receive a "kick" during the explosion. Assuming that ~ 1 Mq of material 
with a kinetic energy of ^ 10^"'^ erg is below the shock when the explosion starts, a ~ 2% asymmetry in the bulk 
motion of this material would result in a kick velocity of ~ 0.02(10^^^ erg/1 Mq)^ ~ 140 km s~^ for the protoneutron 
star. This could explain the observed velocities for a large fraction of pulsars p^q . 

While the angular dependence of the heating rate has some interesting dynamic effects as discussed above, it 
only introduces minor perturbation on the position of the gain radius. To very good approximation, one may use 
the angular average of the heating rate [obtained effectively by dropping the terms involving cos 9 in Eq. 141|l ] in 
determining the gain radius. The angle-averaged heating rate for Y^, = 0.5 is shown as a function of r (solid curve) 
in Fig.|^ In the same figure, we also show the cooling rate qc{B = 10^^ G) as a function of r using T{r) in Eq. 1)16(1 
and Ye — 0.5 for S* = 10 (dashed curve) and S* = 20 (dot-dashed curve), respectively. By comparing Figs. El and |51 it 
can be seen that the gain radius decreases significantly from 137 km for i? = to 99 km for B — 10^^ G in the case 
of S* = 10 but essentially remains at 68 km in the case of 5* = 20. This is because the magnetic field greatly reduces 
rje (see Fig. 0Jd), and hence, the cooling rate (see Fig. |H^) for S — 10. But for S = 20, rje is already small for B — 
(see Fig. 0]d) and reduction of rje by the magnetic field does not change the cooling rate significantly (see Fig. |H^) . 
Numerical models 0,01 show that the material below the stalled shock initially has S ^ 10. Taking T{r) in Eq. H16|l . 
Ve = 0.5, and S = 10, we calculate the gain radius as a function of B and show the results in Fig.^| It can be seen 
that magnetic fields of B ~ 10^^ G or larger significantly decrease the gain radius, thereby enhancing the net heating 
below the stalled shock. Consequently, the shock may be revived more efficiently (i.e., within a shorter time) to make 
an explosion. 



The heating and cooling rates in a uniform field discussed in Sec. IIIII can be adapted to the case of a dipole field 
in a straightforward manner. The strength of the magnetic field to be used in the expressions for x [Eq. I|22|l ] and 
Ee{n,pez) [Eq. (|^ ] is now 



B. Dipole Field 



As a second example, we consider the magnetic field of a dipole in the z-direction: 




(43) 




(44) 
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FIG. 10: The gain radius Rg as a function of magnetic field strength B for a uniform magnetic field (Ye ~ 0.5 and S = 10). 
The gain radius for B ~ 10^^ G approaches that for B = 0. 




FIG. 11: The heating rate per nucleon [dashed {cos 9 — 0) and solid (cosS — 1) curves] and the cooling rate per nucleon 
[dot-dot-dashed {cos6 = 0) and dot-dashed {cos 6 = 1) curves] in a dipole magnetic field as functions of radius r (Ye ~ 0.5 and 
S = 10). The dipole field has strength B{r,e) = Bo{R,./rfVTT3co^ with Bo = 5 x lO" G. 



In addition, the integration over the neutrino sohd angle (see Fig. to obtain the heating rate is changed to 

/ cose^df!, = 2n{R,/rf^J^^L=. (45) 
J Vl + 3cos-^6' 



For illustration, we take i?o = 5 x 10^^ G and show the heating rate as a function of r for cos 9 — (dashed curve) 
and cos6 — 1 (solid curve), respectively, in Fig. ^] In the same figure, we also show the corresponding cooling rate 
using T{r) in Eq. lfTC|l . Yg = 0.5, and 5' = 10 [dot-dot-dashed (cos6' = 0) and dot-dashed (cos0 = 1) curves]. It can 
be seen that the gain radius differs significantly for cos 9 = and 1. For a close examination, we show the gain radius 
Rg as a function of cos 9 (solid curve) in Fig. ^1 Compared with Rg — 137 km for i? = (dashed curve), Rg is 
substantially reduced to ~ 80 km close to the north and south poles of the magnetic field (| cos9\ ^ 1). This reduction 
in Rg is somewhat larger than that for the uniform magnetic field discussed above. This is because the strength of 
the dipole field at r 100 km in the polar directions is somewhat larger than the strength of i? = 10^^ G taken for 
the uniform field. By comparison, Rg only decreases slightly from 137 km for i? = to 131 km near the equator of 
the dipole field (| cos9\ <C 1). This is because at a given r, the strength of a dipole field for | cos9\ <C 1 is weaker than 
that for I cos^l 1 by a factor of ~ 2 [see Eq. (031)]. 

The gain radius for [ cos9\ ~ 1 is significantly smaller than that for [ cos6'| ^ 1 in the case of a dipole field. This 
indicates that the explosion is very likely to occur first in the polar directions of the magnetic field. Note that the 
ejecta in the north and south poles tend to kick the protoneutron star in opposite directions. However, there would 
be a net kick as the heating rate differs for cos 6* = — 1 and 1. This can be seen from Fig. ^1 which shows a small 



14 



140 



130 
120 

B 

A 110 
100 
90 



""-1 -0.5 0.5 1 

cosG 

FIG. 12: The gain radius Rg as a function of cos 61 (solid curve) for the dipole magnetic field in Fig. 1111 Compared with 
Rg = 137 km for B = (dashed curve), Rg is substantially reduced to ~ 80 km close to the north and south poles of the 
magnetic field (] cos9\ ~ 1). 

difference in the gain radius for 008 6* = — 1 (Rg = 80 km) and cos6' — 1 {Rg = 83 km). Thus, we expect that the 
possible kick received by the protoneutron star is similar for both the uniform and the dipole fields discussed here. 

V. CONCLUSIONS 

We have calculated the rates of heating and cooling due to the neutrino processes in Eqs. and ^ in strong 
magnetic fields. We find that for B ~ 10^^ G, the main effect of magnetic fields is to change the equations of state 
through the phase space of e~ and e''^, which differs from the classical case due to quantization of the motion of 
e~ and e"*" perpendicular to the magnetic field. As a result, the cooling rate can be greatly reduced by magnetic 
fields of ^ 10^^ G for typical conditions {S ~ 10) below the stalled shock and a nonuniform protoneutron star 
magnetic field (e.g., a dipole field) can introduce a large angular dependence of the cooling rate. In addition, strong 
magnetic fields always lead to an angle-dependent heating rate by polarizing the spin of n and p. The decrease in 
the cooling rate for magnetic fields of ~ 10"'^^ G decreases the gain radius and increases the net heating below the 
stalled shock. We conclude that if magnetic fields of i? ~ 10^^ G exist within ~ 100 km of the protoneutron star, the 
shock can be revived more efficiently (i.e., within a shorter time) to make an explosion. In addition, the anisotropy 
in the heating rate induced by strong magnetic fields and that in the cooling rate induced by strong nonuniform 
(e.g., dipole- like) magnetic fields may lead to significant asymmetry in the bulk motion of the material below the 
stalled shock, eventually producing an asymmetric supernova explosion. We speculate that this may be one of the 
mechanisms for producing the observed pulsar kick velocities. Obviously, a full treatment of magnetic fields during 
the formation of a protoneutron star and during the supernova process in general greatly increases the complexity 
of an already difficult problem. Nevertheless, we hope that the interesting effects of strong magnetic fields discussed 
here would help to motivate the eventual inclusion of magnetic fields in supernova models. 
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APPENDIX A: EQUATIONS OF STATE 

The relevant equations of state concern the net electron number density (ng- —n^+) [Eq. H17() ] and the total entropy 
per nucleon S [Eq. H18|) ]. The photon contribution to S is: 
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For the conditions in the region of interest, nucleons are nondegenerate and nonrelativistic. So their contribution Sn 
to S is: 



S 



N 



Y„ In 



= 11.5 + In 



3/2' 



2mN f mNT\ 



3/2' 



/ T 1^10^ gem 

I MeV 



Y i^ lOSgcni-3 ^ 



-Y,\nY,-{l-Y,)\r,{\-Y,), 



(A2) 



where we have used y„ = 1 — Yg and Yp = Y^. The above expressions for and are vahd for both the case of no 
magnetic field and the case of magnetic fields of S 10^^ G considered here. 

In the absence of magnetic field, the general expressions for the number densities n^i , energy densities ee± , and 
pressure P^± of e"*" and e~ are: 



■7r2 , 
^ -10 



2^3 



exp[(i;e±/T)±77e] + l 

> rr.2 



exp 



+ (me/T)2 ± ,y. 



■ dx, 



+ 1 



(A3) 



and 



TT^ Jo exp[(i;,±/T)±77e] + l 



■4 /"oo 



exp 



■ dx, 



+ 1 



avr^io exp[(£;,±/r)±,7,] + l 
3^ 



exp 



■ dx. 



(A4) 



(A5) 



In Eqs. (EU-lEni), the upper sign is for and the lower sign for e . 

In magnetic fields, the energy levels and the phase space of e+ and e~ are changed [see Eqs. (|25ll and H27() ]. The 
corresponding expressions for n^± , Cgi , and P^± are: 



^^0^"Jo exp[(i?e±/r)±77e] + l 



n=0 



1 



dx. 



exp 1^x2 + [{ml + 2neB)/T^] ± Ty^j + 1 



(A6) 



/"^ E,± 



eBT^ g ^„ r V.2^^(,,2^,.,^)/^ 2j ^^^^ 



27r2 

n=0 



exp |v^a;2 + {[ml + 2neB)/T2] ± ?7e j 



and 



2n 



^-^t,.r /N^^^yirnl^^r.eB)/ r^\ a.. (A8) 

„=o -^0 exp |v/a;2 + [(m2 + 2neB)/T^] ± rj^j + 1 
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In Eqs. (jA6|) - (|A8|) . the upper sign is for e"*" and the lower sign for e . 

The contributions Sf,± from and e~ to S can be obtained in terms of n^± , £g± , and P^± as 



More specifically, 

where we have used Uf.- — nf.+ — Yep/rriN. 



P.- + P. 



T{p/mN) 



(A9) 



(AlO) 



APPENDIX B: NEUTRINO PROCESSES IN STRONG MAGNETIC FIELDS 



A number of studies on neutrino processes in strong rnagnetic fields exist in the literature. The forward and reverse 
processes in Eq. have been studied in Refs. and [20| assuming that e~ and p are in the ground (n = 0) Landau 
levels. All the four processes in Eqs. and (jSjl have been studied in Ref. assuming that e~ and only occupy 
the ground Landau levels. The forward processes in Eqs. ^ and l(2Jl have been studied in Refs. 33] and j23i] assuming 
that the magnetic field only affects the phase space of e~ and e"*". Parity violation in the forward process in Eq. 
has been studied in Ref. [2j. The cross section of the forward process in Eq. has been calculated in Ref. |23 
using an approach similar to ours. In this appendix, we treat the forward and reverse processes in Eqs. and (O 
in magnetic fields of i? '--^ 10^^ G. Such fields are strong enough to change the motion of e~ and e+ but do not affect 
the description of weak interation (see Sec. IIII|I . 

Ignoring higher order corrections, we take the effective four-fermion Lagrangian of weak interaction to be 



Gp 

71 



cos 6lr7 A^IL" + h.c., 



(Bl) 



where h.c. means the Hermitian conjugation of the first term, 
classical V ~ A form 



In Eq. (|B1|1 . the leptonic charged current has the 



(B2) 



and the nuclconic current N°' is 



(B3) 



The form factors / and g in Eq. IIB3p are taken as constant. In the calculation below, we shall use the Dirac-Pauli 
representation and take the magnetic field B to be in the positive z-direction. All the terms of order and higher 
are ignored in the calculation. 

The wavefunction of a left-handed neutrino with momentum = i?,y(sin G^^, 0, cos 0^) is 



oi(p„-x-B„t) 

\/2L3/2 



/ sii 



sm 



el 



oi(Pi,-x-iS„t) 

- U" 

L3/2 



(B4) 



V cos%^ y 



where L is the linear size of the normalization volume. The above wavefunction also applies to a right-handed 
antineutrino with the same momentum. The wavefunction of a non-rclativistic nucleon is 



gi(k„(p) ■x-m„(p)t) 



V ) 



gj(k„(p)-x-m„(p)t) 



L3/2 



n(p) 



(B5) 



where s„(p) = ±1 denotes the spin state and k„(p) is the nucleon momentum. 
In cylindrical coordinates 0, z), the wavefunction of an electron is 
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where n is the quantum number of the Landau level, is the quantum number of the gyromotion center, A = Xj \J eB 
is a characteristic length scale defined by the strength of the magnetic field, and Se = ±1 corresponds to the two spin 
states when the electron is at rest. The spinor V'n.n^.pe^.se ™ 



/(me + ^e)e-''^/„-i,„,(e72A2)\ 




V 



(B7) 



/„,„„(eV2A2) y 



and 



/ 

(TOe+i;e)/„.,v(eV2A2) 

_^v|^e-^*/„_i,„,(CV2A2) 

V -Pe./n.„.(eV2A^) y 



(B8) 



The special function /n,n,. ('y) in the above equations is defined in Ref. |26|. and can be written in terms of the 
generalized Laguerre polynomial L"^"''(77) as 




(B9) 



The electron wavefunction discussed above is the same as given in Ref. up to a phase factor [Note that there are 
a few typos in that reference: Eq. (45) should read ■4'n,i = (— i)"(2"+'+^7rZ!n!)~ 2 A""^ x • • • and Eq. (46) should read 
i^n,i = i"(2"~'+^7ri!n!)~^A~-^(p/A)"~' x • • • ]. The wavefunction of a positron is 



where 



and 



( 



V27rA2L 

/i^e-'^/„_i,„,(eV2A2)\ 

_ -Pe./n,„.(eV2A2) 



V (me + Se)/„,„.(eV2A2) / 



(BIO) 



(Bll) 



1 



V 



*^/„,„.(e/2A2) 
-(me + £;e)e-'^/„-i,„,(eV2A2) 




(B12) 



All the wavefunctions arc normalized to have one particle in a volume of L^. 
The scattering matrix of + ri ^ + p is 



Gf cos 6'(7 

Gfcos9c 1 



\/2 \/2^L5 



<,.,7a(/ - 575)V'L,s„C.,P..,.e7"(l - 75)V';„ d^x 



In Eq. imi, 



971: 



ede / e™---e-("-"'-)^C7£^_,^7a(/ - 9l5)K^,sUln.,p..,s.ril - 75)C/p^ 



(B13) 



(B14) 



where the subscript "_L" denotes a vector in the xy-plane and = (k„ + Pi, — kp)^. Similarly, the scattering matrix 
oi + n —^ + p is. 



Gp cos 6c 
V2 

GpcoaOc 1 



V2 V2^L 



<„.,7a(/ - 375)^L,.„C!n.,P„,..7"(l - 15^1 d^a. 

27r<5(^;^ - i^e - A)27r<5(fcp^ +p^, - - Pez)Tl' . 



(B15) 
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In Eq. ||BT5)| 



2tv 



,7"(l-75)t/p"„ 



(B16) 



where w'^ — (k„ — kp — Pu)j_- The scattering matrices of + p ^ n + ly^ and f'e + p — s- n + e"*" are the Hermitian 
conjugate of those in Eqs. (jB13|) and (IB15|) . respectively. 
Based on formula 8.411.1 in Ref. 27], we obtain 



^iw±p cos{(j)—(j)Q)—i[n—nj^)(f> 



(B17) 



where (jjy, is the azimuthal angle of w_l, and Jniv) is the Bessel function of order n. Using this and formula 7.421.4 
in Ref. 



/•oo / 2 \ 



one can prove that 



(B18) 



CdC 



= i"+"-e-*("-"")'^''47rA2/„,„^ (2A2w2^). 
Using Eq. (|B19p . we obtain \dJl\^ of the forward and reverse processes in Eq. ^ as 
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(B19) 



(B20a) 

(B20b) 
(B20c) 



(B20d) 



where Vez — Pez/ Ee is the longitudinal velocity of the electron. The corresponding expressions of , which apply 
to the forward and reverse processes in Eq. (0), can be obtained from the above expressions of |9Jt| by changing the 
signs of the terms proportional to (/^ — g^) and replacing wj^ with w^. 
The cross section of i^e + ^ + p is 



Ldpe 



1 Gl cos2 Oc 



27r L-^L 



3r-3 



X27r(5(i;e - El, - A)2TTS{kpz + Pez - Kiz - Pyz) 



(B21) 



In this appendix only, the symbol T deontes the duration of a process such as I'e + — > e +p. Evaluation of au^n{B) 
can be simplified using the summation rule [2^ for the special function /„^„^(2A2w^) in |9Jl|2, 



(B22) 
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where 5n,n' is the Kronecker delta function. A difficulty in evaluation of a^^n{B) is that the integrand in Eq. HB21|I 
is independent of (kp)_L in the infinite nucleon mass limit and the integral diverges. Another difficulty is that there 
is a remaining factor of L^^. These difficulties arise because e~ does not have definite transverse canonical momenta 
and we drop all the terms of order mj} and higher. However, e~ has definite transverse kinetic momentum squared 



(7re)i = (pe 
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2A2 



Pey 



2A2 



= (2?i + l)eB. 



The limit on x and y then corresponds to a limit of L/2X^ on p^x and Pey Thus, we take 
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Using Eqs. P20|l . (|B2T 



(2^)2 

221). and (IB24II. we obtain 
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(B26) 
(B27) 



and Xn is the net polarization of the neutron in the initial state. In Eq. (|B26|I . gn denotes the degeneracy of Landau 
level n for e~. The cross section of + P — > + can be obtained in a similar way as 
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As e~ does not have definite velocity, we calculate the volume reaction rate instead of the cross section for e~ +p 
+ n. To illustrate the dependence on the direction of the outgoing i/g, we first calculate the differential volume 
reaction rate 
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1 G^, cos^ ec 



(27r)3 



In Eq. l\Bm 
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x2'kS{E^ + A - Ee)2'K6{puz + Kz - Pez - ^P^ \^X'^L^° 1U 



El ^ — T 1 ^ — > SttA^ ^ — ^ 



(B29) 



(B30) 



which represents the average over all possible initial e states with a given energy E^ and a given Landau level 
quantum number n. Using Eqs. (|B20() . HB24|I (with kpx(y) replaced with knx{y)), HB29|I . and 1B30|I . we obtain 



dre-p(B) 

dcos 



.(0) 

■ ejV 

2 



^5,1,0- 



l + 2xp4t^cose„ 



.(0) 

eN 



■p + 3g^ 
i^cose +2.il±^ 



/2 + 3g2 



(B31) 
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where 



2n 



is the volume reaction rate without magnetic field. Integrating (|B31|I over cosOi,, we have 



T , . O U +9)9 
l+dn,0 iXv « , o„2 



72 + 3.g2 

The volume reaction rate of e"*" + n — + p can be obtained in a similar way as 



dre+„(5) ^ r, 

dCOS 



(0) 



l + 2x 



.(0) 



and 



r,+„(i?) = r(°^ 



(/ + g).9 



COS 81, 



COS 81/ + 2x 



1 + (5„,o 2Xr: 



/2 + 



(B32) 
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